Measurement of the generalized form factors 
near threshold via — > mr + at high Q 2 
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We report the first extraction of the pion-nucleon multipoles near the production threshold for 
the n-K + channel at relatively high momentum transfer (Q 2 up to 4.2 GeV 2 ). The dominance of 
the s-wave transverse multipole (_Eo+), expected in this region, allowed us to access the generalized 
form factor G\ within the light-cone sum rule (LCSR) framework as well as the axial form factor 
Ga- The data analyzed in this work were collected by the nearly 4-7T CEBAF Large Acceptance 
Spectrometer (CLAS) using a 5.754 GeV electron beam on a proton target. The differential cross 
section and the -k — iV-multipole Eo+/Gd were measured using two different methods, the LCSR 
and a direct multipole fit. The results from the two methods are found to be consistent and almost 
Q 2 independent. 

PACS numbers: 



I. INTRODUCTION 

Pion threshold photo- and electroproduction has a long 
history with continuous interest from both experimen- 
tal and theoretical sides. These studies are of interest 
because the vanishing pion mass approximation in chi- 
ral symmetry, supplemented by current algebra, allows 
exact predictions to be made for the threshold cross sec- 
tions, so-called Low-Energy Theorems (LET) (HUH. As 
a prominent example, the LET establish a connection be- 
tween charged pion electroproduction and the axial form 
factor of the nucleon. In the real world, the finite pion 
mass cannot be ignored (m^/m^r ~ 1/7). The study of 
finite pion mass corrections to LET was a topical field 
in high energy physics in the late sixties and early sev- 
enties before the discovery of Bjorken scaling in Deep 
Inelastic Scattering (DIS) and the advent of Quantum 
Chromodynamics (QCD) [||. (A monograph on pion- 
electroproduction [5[ addresses many of these develop- 
ments.) 

In the eighties and nineties, a renewed interest in 
threshold pion production was triggered by the exten- 
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sive data that became available on_7p — > 7r°p and 
7 *p -> n°p at Q 2 = 0.04-0.1 GeV 2 || (q = p e -p e ,,Q 2 = 
—q 2 )- At the same time, the advent of CHiral Perturba- 
tion Theory (CHPT) has allowed the systematic expan- 
sion of physical low-energy observables in powers of the 
pion mass and momentum. The new insight brought by 
CHPT calculations is that certain loop diagrams produce 
non-analytic contributions to scattering amplitudes that 
are lost in the naive expansion in the pion mass Q @ ■ The 
expansion at small photon virtualitics Q 2 has to be done 
with care as the limits m w — > and Q 2 — > do not com- 
mute in general (Io| . The LET predictions that include 
CHPT corrections seem to be in good agreement with 
experimental data on pion photoproduction [ll|. Ex- 
perimental results on the s-wave electroproduction cross 
section for Q 2 ~ 0.1 GeV 2 are also consistent with 
CHPT calculations when chiral loops are taken into ac- 

30. 



count 



We report the extraction of the multipole E Q+ near 
pion threshold in the charged single pion electroproduc- 
tion channel (ep — > e'nir + ) with a nearly 6 GeV electron 
beam incident on a proton target. This experimental 
data set allowed us to study near-threshold pion produc- 
tion at photon virtualities Q 2 up to ~ 4.2 GeV 2 . This 
experiment is a major step forward and requires very 
good energy resolution in order to approach the pion 
production threshold, where the p-wave contribution of 
the Mi + multipole is suppressed. 
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II. LIGHT CONE SUM RULE MODEL 

In the traditional derivation of LET using the Par- 
tially Conserved Axial Current (PCAC) approximation 
and current algebra, Q 2 is not assumed to be small and 
the expansion in powers of the pion mass involves two pa- 
rameters: TOtt /tojv and m^Q 2 jvr? N [J]|9(. At high Q 2 , the 
second parameter can still be kept small but, in this case, 
the pion is not soft in the target rest frame, even though 
at threshold it is soft in the ir — N final state center-of- 
mass frame. For the threshold kinematics, this affects in 
particular the contribution of pion emission from the ini- 
tial state p3< The LET is formally valid (modulo CHPT 
loop corrections [13]) for momentum transfers as large as 



Q 



l N- 



However, no dedicated experimental study 



of threshold pion production in the Q 2 ~ 1 GeV 2 region 
has been carried out so far. For m n Q 2 /m 3 N — 0(1), the 
LET breaks down: the initial state pion radiation oc- 
curs at time scales of the order 1/mjv rather than l/m n , 
requiring additional contributions of hadronic intermedi- 
ate states other than the nucleon. Finally, at very large 
momentum transfers (Q 2 3> 1 GeV 2 ), one can factorize 
hard-scale contributions as coefficient functions in front 
of soft contributions involving small momenta, allowing 
the use of current algebra (or CHPT) for the latter, but 
not for the amplitude as a whole [14j . 

For asymptotically large Q 2 , the standard perturbative 
QCD (pQCD) collinear factorization technique [Hj][l6| 
becomes applicable, and the helicity-conserving Eq+ mul- 
tipoles can be calculated (at least for = 0) in terms of 
chirally rotated nucleon distribution amplitudes. Resolv- 
ing the onset of the pQCD regime at very large momenta 
is difficult because of the competition between the fac- 
torized contribution, proportional to (a s (Q)/2ir) 2 (which 
runs with Q ), and the non-perturbative soft contribu- 
tions. The latter are nominally suppressed by extra pow- 
ers of Q 2 but are susceptible to end-point contributions 
that cause them to contribute even at very high Q 2 . 

The light-cone sum rule (LCSR) approach has ad- 
dressed this problem, and a realistic QCD-motivated 
model for the Q 2 dependence of both the transverse Eq + 
and the longitudinal Lq + s-wave multipoles has been de- 
veloped (nf for the Q 2 ~ 1 - 10 GeV 2 near-threshold 
region that can now be accessed by the presented exper- 
imental data. A technique was developed [l8[ to calcu- 
late baryon form factors for moderately large Q 2 using 
LCSR [lj|[2(|. The same technique PH is applied to 
pion electroproduction. This approach is attractive be- 
cause in LCSR, soft contributions to the form factors are 
calculated in terms of the same nucleon distribution am- 
plitudes that enter the pQCD calculation without dou- 
ble counting. Thus, the LCSR provide the most direct 
relation of the hadron form factors and nucleon distribu- 
tion amplitudes that is currently available, without using 
other nonperturbative parameters. 



The relevant generalized form factors were estimated 
in the LCSR approach [2l| for the range of momentum 
transfers Q 2 ~ 5 — 10 GeV 2 . For this work, the sum 
rules have been re-derived in [2lj . taking into account 
the semi-disconnected pion-nucleon contributions in the 
intermediate state. The applicability of the sum rules 
can be extended to the lower Q 2 region, and the LET 
are indeed reproduced at Q 2 ~ 1 GeV 2 to the required 
accuracy O^m^)- The results presented here essentially 
interpolate between the large Q 2 limit considered in [2l[ 
and the standard LET predictions at low momentum 
transfers. Two generalized form factors are introduced 
that describe the Q 2 -dependence of the s-wave multipoles 
of pion electroproduction at the threshold kinematics. 
In a simplified approach, the LCSR correlation function 
for the electroproduction close to threshold shows 
dominance of s-wave form-factor-like contributions. 



III. THE GENERALIZED FORM FACTORS Gi 
AND G 2 FROM LCSR 

A. Differential cross sections and form factors 

In the one-photon-exchange approximation, the single 
pion electroproduction cross section factorizes as 
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d 2 a v 
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= ot + €(Jl + cott cos 20* 



+ y/2e(l + e)a LT cos 4>l . 

The parameter e represents the virtual photon polariza- 
tion and r„ is the flux of virtual photons. Ei and Et 
are energies of the initial and scattered electrons respec- 
tively. The angle (f>* is the azimuthal rotation of the mr + 
plane with respect to the electron scattering plane (e, e'), 
v (=Ei — Ef) is the energy transfer of electron, 9 e is 
the polar angle of the scattered electron in the Lab sys- 
tem, f2* is the solid angle of pion in the center-of-mass 
frame, and W is the invariant mass. In the absence of 
a transverse polarization of the target nucleon, the cross 
section does not depend on <f> e . For an electron beam and 
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a proton target, the center-of-mass differential cross sec- 
tion d 2 a u depends on the virtual photon polarization (e) 
through four structure functions : ot + £<Jl and the in- 
terference terms <jtt and olt- Four structure functions 
are determined by a fit to the 0*-dependent differential 
cross section. The partial wave decomposition by Legen- 
dre polynomials of the structure functions in the limit of 
angular momenta I < 2 is given by [l7| 

n 

a T + ea L = Y J Dj +L Pi(cosei) , 

1=0 
n-2 

ott = sin 2 6* Dj T Pi(co&6l) , and 

1=0 
n-l 

olt = sine; Df T Pi{cos6l) , 

1=0 



where the coefficients {D^ +L , D^ +L , D%+ L , Dq T , Dq T 
and D^ T ) depend on seven complex multipoles since the 
Legendre coefficients are directly related to multipole de- 
composition. The quantity 9* is the tt + polar angle in the 
center-of-mass frame. In the LCSR approach in the pion 
threshold region with vanishing pion mass, the Legendre 
coefficients can be described in terms of the generalized 
form factors by 
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where G^ and G E are the magnetic and electric Sachs 
form factors of the neutron (due to pion emission off the 
initial proton), = \/2 is the isospin factor, f n = 93 
MeV is the pion decay constant, and qa = 1.267 is the 
axial coupling. For charged pion production additional 

contributions, G™ arise from the chiral rotation of the 
electromagnetic current. These are not present for neu- 
tral pion production. The Legendre moment Dtt is zero 
since d- waves are absent. Paramctrizations developed 
previously [22| for the electric and magnetic neutron form 
factors (G^f and G E , respectively) are used in this anal- 
ysis. 



B. dependence on Gi and G2 



As described above, experimental data for G E at high 
momentum transfers are lacking. The quality of predic- 
tions for G E at high momentum transfers also remains 
poor. Thus, the generalized form factor will be extracted 
here under the assumption of m w ~ with and without 
taking G E into account. At first, if we take G E ~ 0, 
Eq. ^ can be re-written as 
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Here fcj, kf are the center-of-mass momenta in the initial 
and final states, respectively [l7[. The low-energy theo- 
rems (LET) relate the s-wave multipoles, or equivalently 
the form factors G\ and G2 at the pion threshold, to the 
nucleon electromagnetic and axial form factors for the 
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is negligible. Therefore, 
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since G% ~ in Eq. ©, G^ T 

only two terms survive from Eq. since in the mr 

(due to 

absence of d- waves). 
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In this case, we can extract the Eq7 



amplitude by using 
its relation to the form factor (G" 7r+ ) [TtJ both normal- 
ized by the dipole form factor (Go = 1/(1 + Q 2 /no) 2 , 
Mo = 0.71): 
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Alternatively, we could take a non-zero value of G^ into 
account, but there are no constraints on the imaginary 
parts of Gi and G2, since only electron- helicity indepen- 
dent data are available. The real parts of Gi and G2 can 
still be determined from D^ +L and Dq t , but to solve for 
the imaginary parts of G\ and G2, further assumptions 
are required. 



IV. EXPERIMENT 

The measurement was carried out with the CE- 
BAF Large Acceptance Spectrometer (CLAS) [23|]. A 
schematic view of CLAS is shown in Fig. [T] CLAS uti- 
lizes a magnetic field generated by six flat superconduct- 
ing coils (main torus) in an azimuthally symmetric ar- 
rangement. The coils generate an approximately toroidal 
field distribution around the beam axis. The six sectors 
of the magnet are independently instrumented with 34 
layers of drift cells for particle tracking, plastic scintilla- 
tion counters for time-of-flight (TOF) measurements, gas 
threshold Cherenkov counters (CC) for electron/pion sep- 
aration and triggering purposes, and a scintillator-lead 



sampling array (electromagnetic calorimeter or EC) for 
photon and neutron detection, as well as triggering. To 
aid in electron/pion separation, the EC is segmented into 
an inner part facing the target and an outer part further 
away from the target. 

CLAS covers on average 80% of the full 4ir solid angle 
for the detection of charged particles in the laboratory 
frame. Azimuthal angle acceptance is maximum at large 
polar angles and decreases at forward angles. Polar angle 
coverage ranges from about 8° to 140° for the detection 
of 7T+ . Electrons are detected in the CC and EC for polar 
angles from 15° to 55°, with this range being somewhat 
dependent on the momentum of the scattered electron 
and the magnetic field strength. 

The target is surrounded by a small toroidal magnet 
(mini-torus) with non-superconducting coils. This mag- 
net is used to shield the drift chambers closest to the 
target from the intense low-energy electron background 
resulting from M0ller electron-scattering processes. In 
the current experiment, only two charged particles need 
to be detected, the scattered electron and the produced 
7r + , while the full final state is reconstructed using four- 
momentum conservation constraints. The continuous 
electron beam provided by CEBAF is well suited for mea- 
surements involving two or more final state particles in 
coincidence, leading to very small accidental coincidence 
contributions of < 10~ 3 for the instantaneous luminosity 
of 10 34 cm _2 sec _1 used in this measurement. 

The specific experimental data set used for this analysis 
was collected from October 2001 to January 2002, during 
the CLAS el-6 run period. The incident beam had an 
average intensity of 7 nA and an energy of 5.754 GeV. 
The liquid hydrogen target was 5 cm long and located 
4 cm upstream from CLAS center. The main torus mag- 
net was set at 90% of its maximum field. Empty-target 
runs were performed to measure contributions from the 
target cell windows. We compared our optimized beam 
energy of 5.754 GeV with the energy calibrated by Hall- 
A, which is based on concurrent high-resolution elastic 
electron-proton scattering measurements. Both beam en- 
ergies agree within less than 6 MeV. 

Raw data were subjected to the calibration and recon- 
struction procedures that are part of the standard CLAS 
data-analysis chain. The reaction studied in this paper 
contributed to only a small fraction of the total event 
sample, and a more stringent event selection ("skim- 
ming") was applied to select events with one electron 
candidate and only one positively charged track. These 
events were subject to further selection criteria described 
in the following sections. The kinematic range and bin 
size were optimized according to the available statistics 
in the covered kinematic range. Table U reports the kine- 
matic range, bin size, and number of bins for the rele- 
vant variables. Acceptance and radiative corrections were 
specifically calculated for each bin given in Table HI The 
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Large-angle Calorimeter 
Drift Chambers I Electromagnetic Calorimeter 




FIG. 1: (Color online) Schematics of the CLAS detector sys- 
tem. The top panel shows a horizontal cut through sectors 1 
(upper hemisphere) and 4 (lower hemisphere) along the beam 
line. The beam enters from the left into CLAS. A GEANT- 
simulated event is shown with an electron bending towards 
the beam line and a positive particle in the opposite sector 
bending away from the beam. The bottom panel shows a cut 
perpendicular to the beam line through the center of CLAS. 



cross section is calculated by multiplying by the radiative 
correction factors and dividing by the acceptance correc- 
tion factors. The 0* is the polar angle of the detected 
positive pion in the center-of-mass frame. 



TABLE I: The ranges of kinematical bins used in this analysis. 



Variable Number 
of bins 


Range 




Bin size 




W 


3 


1.11 - 1.15 GeV 




20 MeV 




Q 2 


5 


2.12 - 4.16 GeV 2 




variable 




cos 9% 


10 


-1.0 - +1.0 




0.2 






12 


0° ~ 360° 


30° ( 


cos 9% > 


-0.1) 




6 


0° - 360° 


60° ( 


cos 9* < 


-0.1) 



V. DATA ANALYSIS 

A. Particle identification and corrections 

For the particle identification (PID) and kinematic cor- 
rections we applied the standard PID cuts for the near- 
threshold physics regime. The total number of single- 
pion events with W < 1.2 GeV (see Fig. [2J is approx- 
imately 4.55 x 10 4 . Since PID and kinematic correc- 
tions have a strong dependence on event statistics, our 
PID and corrections were investigated before applying the 
W < 1.2 GeV cut to avoid large uncertainties from such 
small statistics. Therefore, most of the correction proce- 
dure for electrons and pions follows the method described 
in the previous analysis (25| with optimized parameters. 



1. Electron identification 

Electrons were tentatively identified in CLAS at the 
trigger level during data acquisition by requiring a mini- 
mum amount of energy deposited in the electromagnetic 
calorimeters (EC) in coincidence with a signal in the 
Cherenkov counters (CC). This tentative identification 
was then improved by applying additional requirements 
in the off-line analysis. Coincident hits between the EC 
and CC were also matched with a track reconstructed in 
a drift chamber (DC) in the appropriate sector of CLAS. 
The direct correlation between the energy deposited in 
the EC and the momentum obtained from the track re- 
construction in the DC was used to remove the residual 
pion contamination. 

About 30% of the total energy deposited in the EC 
is directly measured in the active scintillator material; 
the remainder of the energy is deposited mostly in the 
lead sheets interleaved between the scintillator sheets as 
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cos9 n * 

FIG. 2: (Color online) The kinematical coverage for pions 
used in this analysis in terms of the momentum transfer Q 2 
versus the center-of-mass energy W (top) and 0* versus cos 8^ 
(bottom). 



showering materials. To improve the separation of elec- 
trons and pions, the ratio E tot /p was used, where E tot 
and p are the total energy deposited in the calorimeter 
and the deduced momentum for the particle, respectively. 
This ratio, which is called the sampling fraction (a), is 
nearly momentum independent for the range of electron 
momenta (2.5 to 4 GeV) in this analysis. The sampling 
fraction for electrons (determined by all electrons over the 
full W range) was found to be 0.291 for this experiment, 
a value roughly comparable with the estimate of 0.231 for 
that ratio determined in a Monte Carlo simulation. Fig- 
ure [31 shows the application of the sampling fraction cut 
to experimental (left) and simulated (right) data. The 
solid lines represent the ±3cr sampling fraction cuts for 
the experimental data and the Monte Carlo simulation in 
this analysis. 

Most of the produced pions that passed through the 
EC were minimum ionizing particles. These pions de- 
posited energy both in the inner and the outer stacks of 
the calorimeter in amounts almost independent of their 
kinetic energy and only related to the detector thickness. 
Pions were identified by the precise correlation between 
the minimum total energy deposited in the calorimeter 
Etot and the energy deposited in the inner part of the 
calorimeter Ej n . To avoid continuous triggering on noise, 
a minimum signal threshold was set for the calorimeter. 
For this experiment, the threshold was set at a level of 
172 mV, such that only electrons with momenta greater 
than about 640 MeV were detected. 

A GEANT simulation (GSIM) was used to determine 
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FIG. 3: (Color online) Optimized sampling fraction constant 
of EC versus electron momentum with W < 1.2 GeV for the 
experimental data (left) and Monte Carlo simulation (right). 
The solid lines (left) show the ±3er cuts from all electron fit. 
PI is the sampling fraction from fit. 



the response of the electromagnetic calorimeter as a func- 
tion of electron energy. When an electron hit is close to 
the calorimeter edges, part of the shower leaks outside 
the device; in this case, the energy cannot be fully recon- 
structed from the calorimeter information alone. This 
problem was avoided by selecting only those electrons ly- 
ing inside a fiducial volume within the electromagnetic 
calorimeter that excludes the detector edges. 

Tracking information from the drift chambers was used 
to reconstruct for each event an originating vertex loca- 
tion in the target region. Particle identification was im- 
proved by eliminating events from the analysis that had 
reconstructed vertex positions outside the known volume 
of the target (which included a small target misalign- 
ment from the beam axis). For this experiment, these 
vertex requirements demanded that the reconstructed z- 
vertex position (distance along the beam axis from the 
center of CLAS, with positive values indicating down- 
stream of the center) lie in —80 mm < Z vtx < —8 mm. 
We corrected the x-, y- vertex positions for the beam 
centering on the target, and the z- vertex cut was im- 
posed on the reconstructed vertex locations based on the 
beam axis. Since the beam position was not precisely 
centered on the target, with offsets of X tg t = 0.90 mm 
and Itgt = —3.45 mm, the z-vertex was corrected for 
this small misalignment of the beam position before the 
z-vertex cut was imposed on the reconstructed vertex lo- 
cation. Figure [4] (left) shows the z-vertex distribution 
for sector 3 before and after the vertex correction, and 
the z-vetex cuts that have been applied. Figure 0] (right) 
shows the transverse beam position. We also take an 
empty target contribution into account. 

Coincident hits between the EC and CC were also 
matched with a track characteristic for a negative particle 
that is reconstructed in the drift chambers of the same 
CLAS sector. A lower threshold on the number of pho- 
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toelectrons detected in the photomultiplier tubes of the 
CC for an event provided an additional cut for improv- 
ing electron identification. The number of photoelectrons 
detected in the CC sectors follows a Poisson distribution, 
modified for irregularities in light collection efficiency for 
the individual elements of the array. For this experiment, 
a good electron event was required to have more than 2.5 
photoelectrons detected in CC. 




FIG. 4: (Color online) z-vertex distribution for sector 3 be- 
fore (black solid line) and after (red solid line) vertex cor- 
rection (left). The vertex cuts are illustated by the vertical 
dashed lines. Reconstructed X and Y target positions, show- 
ing an offset of X tg t = 0.90 mm and Ftgt = —3.45 mm, re- 
spectively (right). 



2. Pion identification 

For 7r + identification a track characteristic for positive 
particles has to match with a corresponding hit in the 
TOF system. Pions are then separated from other pos- 
itively charged particles based on their hadron velocity 
Ph = v / c, which is obtained from the difference between 
the vertex start time and the time of flight of the TOF 
counters (SC), and their hadron momentum ph, which is 
determined by tracking the hadron through the magnetic 
field in the drift chambers. To isolate pions from protons, 
a ±2(7 cut on j3h versus pu is applied. 

Figure [5] shows the particle's velocity versus momen- 
tum for positive tracks. A single Gaussian fit function 
was applied to (3h in each momentum bin to choose the 
proper velocity cut for pion identification. The solid red 
lines superimposed on the scatter-plots in Fig. [5] show the 
Ph cuts for both simulation (upper) and data (lower). All 
remaining positrons were considered as pions due to the 
limited momentum resolution, which increases the back- 
ground. However, the missing mass and vertex cuts re- 
duce the background to a few percent. 

We observed that the other positive charged parti- 
cles (protons, kaons) contamination under the missing 



mass peak is negligibly small, even at high pion mo- 
menta. Overall, background contributions under the neu- 
tron missing mass are negligible in our kinematic region; 
hence, no background subtraction was performed. 




4 £ 

[GeV]" 



FIG. 5: (Color online) The particle velocity distribution fJ 
versus particle momentum. The pion velocity sample is fit by 
a function of the form A + exp(B + Cp n ) to generate the fi cut 
for both simulation (upper) and experimental data (lower). 



B. Fiducial cuts 

Due to the complexity of the CLAS geometry and edge 
effects, we define a fiducial volume that is restricted to 
detector regions with nearly full particle acceptance and 
high reconstruction efficiency. We made individual CLAS 
sector dependent geometrical fiducial cuts and applied 
the same cuts to the simulations and to the experimental 
data, where those cuts select areas of uniform detector 
response that can be reproduced by Geant3 based simu- 
lation (GSIM) with CLAS detetor geometric information. 
In order to implement the sector fiducial cuts, the GSIM 
program only requires knowledge of the momentum and 
charge of the particle, since the curvature of the trajec- 
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tory of the particle depends only on the particle's charge 
and momentum, and the strength of the CLAS magnetic 
field. The fiducial cuts have been applied to both elec- 
trons and pious in the same way as in the previous anal- 
ysis [25| using optimized parameters. 



1 . Electron sector fiducial cuts 

Since the toroidal magnetic field bends the electrons in- 
ward, the fiducial cut in polar angle 9 e and azimuthal an- 
gle 4> e depends on the momentum of the electron p e [25[ . 
The azimuthal symmetry of the angular distribution for 
the electrons was considered when selecting these regions. 
Therefore, for fixed 9 e and p e , one expects to find a flat 
distribution in <j) e . Several cut functions to eliminate de- 
pleted regions in momentum p e and angle 9 e were applied 
empirically. Figure [6] shows such a sector centroidal elec- 
tron angle <p e distribution in a sample momentum bin of 
3.02 ± 0.125 GeV. The black solid curve in the left plot 
shows the boundary of the fiducial region for the central 
momentum in that bin. Only events with electrons inside 
the black curve are used in the analysis. In addition, a set 
of 9 e versus p e cuts are used to eliminate the areas with a 
depleted number of events due to bad time-of-flight coun- 
ters, Cherenkov counter photomultiplier tubes, or drift 
chamber wires. 




e e [deg] o e [deg] 



FIG. 6: (Color online) The electron angular distribution in 
a sample momentum bin p e = 3.02 ± 0.125 GeV for sector 4 
(left). The black solid curve indicates the fiducial cut bound- 
ary. The criterion to determine the electron fiducial region in 
terms of <f> e for one momentum bin is the detector efficiency. 
In order to eliminate the depleted regions of the detector we 
chose the flat area by looking through the ^-sliced <f) e distribu- 
tions. The right plots show examples of the <j) e distribution in 
two 8 e bins: 21° and 22° for sector 4. The highlighted areas 
around the center indicate the selected fiducial range. 



2. Pion sector fiducial cuts 

Pion fiducial cuts are designed to exclude regions of 
non-uniform acceptance due to interactions with the 
mini-torus coils, torus cryostat, or the edges of the drift 
chambers. The onset of these regions is not symmetric 
about the sector mid-plane but is sector as well as mo- 
mentum dependent. This cut is a function of the pion 
momentum p v after selecting good electrons. 

The pion momentum was scanned with a 100 MeV bin 
width from 0.3 GeV to 1.7 GeV. The pion angular distri- 
bution was investigated in each momentum bin. The tt + 
fiducial cut functions are parametrized by the pion mo- 
mentum and the angles 9 n and (f> n . The same pion cuts 
are applied to data as well as simulation. To define the 
reasonable detector response region, the two-dimensional 
plot of versus n was sliced along 9^ in 2° bins for 
12° < 9 n < 100°, and each distribution was fit by a 
function which included a trapezoidal shape and a con- 
stant. A fit example is shown in Fig. [7] Each momentum 
and 9 n bin in each sector has a unique plateau region. 
The corresponding fit parameters are functions of sector, 
Ptt, and 6„. The correlation between and 9^ is de- 
scribed by an exponential and a third order polynomial 
function. 




FIG. 7: An example of a (f>^ (deg) fit with a trapezoidal 
function. The flat region was used to select the sector fiducial 
cut for the kinematical bin with 0.35 GeV < < 0.45 GeV 
and 23° < 9 n < 25° bin for sector 4. 



C. Kinematic corrections 

The kinematic corrections from the previous analy- 
sis [25| were applied in this analysis. The corrections 
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are based on the measured angles and momenta of the 
detected particles (electrons and pions). Figure |S] (left) 
shows the neutron missing mass for W = 1.11 GeV, 
which is the lowest W bin that we can access in this data 
set. Overall, there is no systematic (^-dependence on 
the inferred neutron mass. The central horizontal dashed 
line in the right-hand plot of Fig. [5] indicates the mean 
value found for the neutron mass in the left-hand plot; 
the upper and lower horizontal dashed lines indicate ±3cr 
of the missing mass distribution, where a is the energy 
resolution as defined by the standard deviation. 



Mean = 0.936 
icma = 0.01 87 
X 2 =2.4003 




5ECT0R=1 





















f 


1 




f 







1 1.1 1.2 
MMx [GeV] 



20 

. [deg] 



FIG. 8: An example of a neutron missing mass fit for sector 
1 by a single Gaussian function in the bin at cos9n ~ 0.1, 
W = 1.11 GeV and Q 2 = 2.44 GeV 2 (left), and the neutron 
missing mass <j> w dependence for sector 1 (right). Dashed lines 
are defined in the text. 



VI. MONTE CARLO SIMULATIONS 

The detection efficiencies for various particles in CLAS 
were estimated from Monte Carlo simulations of the de- 
tector. We used the two different event generator^ AAO- 
RAD (a physics- model-based event generator [23|) and 
GENEV (a phase space event generator [24[). Both 
event generators generate exclusive events including ra- 
diative effects. The AAO-RAD event generator uses the 
MAID2003 model [H|. The GENEV event generator al- 
lows us to generate several exclusive electroproduction re- 
actions, from pion production to the production of vector 
mesons (cj, p°, and 4>) including their decays, as well as 
non-resonant multi-pion production. Cross section tables 
for these processes are used that are based on photopro- 
duction data and extrapolated to the case of electropro- 
duction. 

Several million (AAO-RAD) and ~ 50 million 
(GENEV) single-pion events were generated near the 
pion threshold region between 1.1 GeV <W< 1.2 GeV 



and 1.0 GeV 2 <Q 2 < 10.0 GeV 2 with full angular cov- 
erage. Both generated data sets were processed by the 
standard GEANT simulation of the CLAS detector. The 
statistical uncertainties for both simulations are less than 
a percent. Most steps closely followed the previous anal- 
ysis described in Ref. 



25] 



The comparison of the detection efficiency results us- 
ing the two different event generators indicates a 4% 
uncertainty in the simulated efficiency prediction. Scat- 
tered electron events were simulated by AAO-RAD and 
GENEV for exclusive single charged pion electroproduc- 
tion on a liquid hydrogen target in CLAS for an incident 
electron energy of 5.754 GeV and a CLAS torus current 
of 3375 A. Figure [5] shows examples of results from the 
simulations of W, Q 2 , and angular distributions. 
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FIG. 9: (Color online) Examples of Monte Carlo simulation 
results: Q 2 versus W and the angular distribution <p% versus 
cos 6% for all W (left column), and yield versus Q 2 and the 
angles 0* versus cos (9* for W < 1.2 GeV (right column) as 
a snap shot after GSIM reconstruction but still without the 
GPP step. The two vertical solid lines on the top-left plot 
show the W region between 1.10 GeV and 1.16 GeV of inter- 
est. 



The GSIM Post-processor (GPP) is used for fine ad- 
justments of the reconstructed GSIM data to better 
match the measured data. GSIM simulates events in 
an ideal detector system and GPP is used to adjust two 
quantities. One is the drift chamber position resolution 
smearing factors that affect the tracking momentum res- 
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olution. The other is the TOF time smearing factor that 
affects the timing resolution. The optimized GPP set- 
tings were taken from Ref. (25j . 

The GPP output has been processed with the event 
reconstruction software. In order to relate experimental 
yields to cross sections, one needs to calculate the accep- 
tance, including the efficiency of the detector and radia- 
tive effects. Both acceptance and radiative corrections 
are processed similarly as in the previous analysis (2f|. 



A. Acceptance correction 



To calculate the acceptance for the ir + n channel in the 
CLAS detector system, one can define the acceptance in 
a given kinematic bin i as 



Acceptance^ = 



Ni 



EEC 



jyGEN 



section and J cr NORAD the un-radiated model cross sec- 
tion with integrated luminosity. N RAD and N^ ORAD 
are event numbers for radiative and non-radiative events, 
respectively. The detailed procedure of radiative correc- 
tions is described in [25| . 

Since there are several models to describe the cross 
sections near the threshold, it is important to verify 
whether ExcluRad gives consistent correction factors in- 
dependent of the physics models. In order to perform 
this stud y, w e used two different models, which are Sato- 
Lee2004 [28[ (dynamic model) and MAID2003 [II] (uni- 
tary isobar model) as ExcluRad inputs. Both cover the 
kinematic region of interest. Figure [TTJ shows examples 
of the radiative-correction comparison between the two 
models, as a function of cos 9* at a fixed W = 1.15 GeV 
and Q 2 = 2.91 GeV 2 but for different 0; = 105°, 135°. 
As Fig. [TTJ shows, the radiative corrections from both 
models are consistent even where their cross sections are 
quite different. 



where N® is the number of events generated and 
jyREC j g tfie number of events reconstructed after im- 
posing all cuts. Figure [10] shows typical examples of (^In- 
dependent acceptances using the AAO-RAD event gen- 
erator. 
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FIG. 10: Examples of the (^.-dependent acceptances for the 
7*p n-K + reaction at W = 1.13 GeV, Q 2 = 2.12 GeV 2 , and 
cos6»; = 0.1 - 0.5. 



0.95 

0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 
0.55 

0.5 



1 -0.5 0.5 

<p„" = 1 G5 Meg] 

A = SLee2004 
- O = MAID2003 



0.5 

COS'$„ 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
. 





= 1 35 [deg] 






= SLee2004 






= MAID2003 




-o a 


@ K is a Q o 




t 


. • • * 

• A * A 
A 


• • 










0.5 0.5 






CO 





FIG. 11: (Color online) Cross sections (solid symbols) 
and radiative correction factors (open symbols) as functions 
of cos (J* in a fixed bin centered at W = 1.15 GeV and 
Q 2 = 2.91 GeV 2 . The two plots are for different values of 
the azimuth, <j>% = 105° and 135°. Each plot shows the RC 
factor on the left j/-axis and differential cross sections on the 
right j/-axis. The (red) triangles (open and closed) are from 
the Sato-Lee2004 model [28| ]. The (blue) circles (open and 
closed) are from the MAID2003 model. 



B. Radiative corrections 



The radiative corrections (RC) were calculated using 
ExcluRad [26]. ExcluRad gives an estimate calculation 
of radiative effects in given kinematic bins. The gen- 
eral definition of the radiative correction factor is the 
ratio of radiative events to events without radiative ef- 
fects in a given kinematic bin. It is defined in a fixed 
W, Q 2 , 61, and ^ bin (j): RC, = {N? AD / f a RAD ) ■ 

(J a NORAD /Nf ORAD ), where RCj is the radiative cor- 
rection for bin (j), f a RAD is the radiated model cross 



VII. RESULTS 
A. Differential cross sections 

In order to extract the differential cross sections, all 
previously described corrections, efficiencies, cuts, photon 
flux and luminosity normalizations have to be applied to 
the data. The azimuthal angle </>* is divided into 60° 
bins for cos#* < — 0.1 (yielding 6 bins) and 30° bins for 
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cos (9* > —0.1 (yielding 12 bins). This binning and a 
minimum acceptance cut (> 0.6%) limit the statistical 
uncertainties. 

Figure [TU shows an example of the 0* -dependent dif- 
ferential cross section (left) and the corresponding fit 
(right) in a W = 1.11 GeV bin that is closest to the 
pion threshold in this analysis, with comparisons to the 
Dubna-Mainz-Taipai (DMT) [13], Sato-Lee 2004 HH, 
and MAID2003 ^ models. The \ 2 is calculated by 
Xo/(N p ts — 3), where \o is not normalized to the degrees 
of freedom, N p t s is the number of non-zero data points 
in the ^-dependent histogram, and 3 is the number of 
free parameters in the fit function Eq. ([7]). The overall 
averaged x 2 is 1.25. 

Table |TT] summarizes the systematic uncertainties in 
this analysis averaged over all accessible kinematic bins 
shown in Figures [1411161 The average total systematic 
uncertainty is about 11.5%. The major sources of the 
systematic uncertainty are the bin effects, the particle 
identification (PID) for clcctron/pion separation and the 
different physics event generator calculation of the ac- 
ceptances. The systematic uncertainty for bin effects 
takes into account bin-size, bin-centering, and different 
bin number of 0* with qudraturc sum. The radiative 
corrections show little systematic dependence on differ- 
ent physics models. 



B. Extraction of structure functions 

The fit of the </>*-dependent cross sections allows us 
to access the polarized structure functions. The fitting 
function has three fit parameters, corresponding to the 
structure functions <jt + e &lt, and <jtt'- 



d 2 a 

dn* 



A + C cos 24>l + D cos < 



(7) 



The relation between the structure functions and fit pa- 
rameters is given by A = ot + e <Jl , C = e utt , and 
D = \/2e(l + e) <tlt [25|. Figure [14] shows the structure 
function of ar + euL as a function of cos 9* in different Q 2 
bins near the pion threshold region. Figures [15] and [TBI 
show the interference terms utt and <jlt, respectively. 
For the three structure functions shown, several features 
are notable. First, the E 0+ multipole in the MAID2003 
model plays a dominant role (compare between red bold 
dash and black dash-dot lines) in both or + and cttt 
in the forward angles, and the lower values of Q 2 show 
larger differences between the results with and without 
Eq + . Second, the experimental results show a good con- 
sistency with the MAID predictions at all Q 2 . Finally, 
most of the ott results particularly at high Q 2 are close 
to zero in the W — 1.11 GeV bin, which is the closest to 
threshold. It is expected that the d-wave contribution is 
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FIG. 12: (Color online) Examples of the differential cross 
sections as function of (fa with 0.6% minimum acceptance cut 
at W = 1.11, 1.13 GeV, Q 2 = 2.12, 2.44, 2.92 GeV 2 , and 
cos#£ = 0.1, —0.3. The error bars of the data include only 
statistical errors and the shaded bars show the systematic un- 
certainties. Results from the DMT model [2tJ are indicated by 
the black line. Also shown are the Sato-Lee 2004 model (red 
dashed-dot line) JH) and MAID2003 (blue dashed line) [2<|. 
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FIG. 13: (Color online) Examples of fits to the differential 
cross sections as function of <j> J at W = 1.11, 1.13 GeV, Q 2 = 
2.12, 2.44 GeV 2 , and cos^ = -0.3, 0.1. The corresponding 
fit is in Eq. Q. 



absent at threshold, although the MAID calculations still 
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TABLE II: Average systematic uncertainties from various 
sources for the differential cross sections from this analysis. 



Source 


Criterion 


Estimated 
contribution 


e~ PID 


sampling fraction cut in EC 
(3cr -»■ 3.5cr) 


4% 


e~ fiducial cut 


width (10% reduced) 


2.2% 


7T+ PID 


P resolution change 
(2cttof — > 2.5<ttof) 


1.3% 


tv~^~ fiducial cut 


width (10% reduced) 


3% 


missing mass 
cut 


neutron missing mass resolution 
(3<tmmx — > 3.5<tmmx) 


1% 


vertex cut 


z-vertex width (5% reduced) 


1% 


acceptance 
correction 


event generator dependence 
AAO_RAD versus GENEV 


4% 


radiative 
correction 


physics model dependence 
Sato-Lee2004 versus MAID2003 


0.5% 


mininum 
acceptance cut 
and bin effect 


aplied 0.6% cut 
and 6,12 (f>% 
bins 


9% 


Total 




11.5% 



indicate substantial d- wave at low Q 2 . 

The LCSR calculations predict ^-independent differ- 
ential cross sections caused by a complete cancellation 
of the contribution to oxt and ott from G\ and G2 
(see Section ITTT|) at particular values of Q 2 . This analysis 
shows the #* dependence for the longitudinal-transverse 
interference term cr iT , as seen in Fig. \TE\ a LT is es- 
pecially strong for large angles and relatively small Q 2 . 
However, this is still quite consistent with the LCSR pre- 
diction if we focus on the highest Q 2 bin, because the <j>* - 
independence is only expected for large Q 2 in LCSR jl7j |. 
The experimental data for cjlt and ott show small de- 
viations from zero over all cos 6** for all three W bins at 
Q 2 = 4.16 GeV 2 , which is possibly caused by the previ- 
ously mentioned cancellation of the G\ and G2 contribu- 
tions. 



C. Extraction of Legendre moments 

The Legendre moments for the components of the dif- 
ferential cross sections were extracted by fitting the cos #* 
distributions of the <tt + &tt, and <tlt structure 
functions with first and second order Legendre polynomi- 
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FIG. 14: (Color online) The structure function (<7t + eo\l) 
as a function of cost9* at W = 1.11 — 1.15 GeV and Q 2 = 
2.12 - 4.16 GeV 2 with model predictions : MAID2003 (red 
bold dash: full multipoles, green bold dash: without So+, an d 
black dash-dot: without E 0+ ) and MAID2007 (blue bold dot). 
The shaded bars show the estimated systematic uncertainties. 



als. As mentioned in Section Ilia (see Eq. [2]), Pi (cos 8*) 



(T+L) n (TT) 



is the I -order Legendre polynomial and D, , D, , 

(LT) 

D) are the Legendre moments for <jt + e<7L, <JtTi and 
CiTj respectively. Each moment can be written as an ex- 
pansion of magnetic (Mi ± ), electric (Ei ± ), and scalar 
(Si ± ) 7rV-multipoles or as a factorization of generalized 
form factors in the light-cone sum-rule framework. The 
expansion is truncated at — 1, because the s-,p-wave 
interference terms, particularly those involving the mul- 
tipole Eq + , dominate near the threshold. 



Figure [T71 shows the Q 2 -dependent Legendre moments 
for ot + e&L, ott and olt at W = 1.11, 1.13, and 
1.15 GeV with different model predictions. Our mea- 
surements cover 1.9 GeV 2 < Q 2 < 4.5 GeV 2 with bin 
center values from 2.12 to 4.16 GeV 2 . Figure [T71 reveals 
that the extraction of Dq T leads to values close to zero 
over all Q 2 , particularly in the lowest W bin, whereas the 
MAID models predict a sizable amplitude. 
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FIG. 15: (Color online) The structure function ott as 
a function of cos#^ at W = 1.11 — 1.15 GeV and Q 2 = 
2.12 - 4.16 GeV 2 with various model calculations. Curves 
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D. Extraction of the Eq+ by the LCSR method 

The extracted Legendre moments can be directly used 
in Eq. ((3]) or Eq. Q depending on the parametrization 
of the electric form factor of the neutron (G%) and other 
kinematic constants. Figure [18] shows the results of the 
Q 2 -dependent Eq+ multipole divided by the dipole form 
factor (Gd) in the W bin nearest to the threshold. The 
plot on the left of Fig. [TS] shows the comparison of Eq + 
with and without taking systematic uncertainties into ac- 
count under the assumption of GJ. = 0. The plot on the 
right shows the comparison of Eq + multiplies extracted 
with G% = or G% ± 0. MAID2007 [H and LCSR 
calculations are also shown. 

The experimental results for Eq + /Gd are about 0.2 — 
0.3 GeV -1 and almost flat as a function of Q 2 . The am- 
plitude is larger than predicted by MAID2007 and sim- 
ilar to (or a bit smaller than) the LCSR calculations, 
although LCSR has a steeper Q 2 dependence, which may 
be caused by the extrapolation in the current LCSR cal- 
culations of the form factors to the chiral limit (m w — > 0). 
The multipole was extracted for two different G% dipole 
parametrizations (3ll . |32[ under the assumption of van- 
ishing pion mass. Figure [19] (left) shows that the com- 



parison of both G% parametrizations leads to negligible 
difference. 

Since CLAS has published the measurement of the 
magnetic neutron form factor at hi gh m omentum trans- 
fers with a level of accuracy of 3% [33], we can directly 
substitute this result into our extraction instead of the 
parametrization. Figure [19] (right) shows the results of 
E 0+ at W = 1.11 GeV for the measured CLAS G n M form 
factor. 



E. Extraction of Eo+ by multipole expansion 



The relations between the invariant amplitudes fi and 
the helicity and multipole amplitudes of the cross sections 
are found in Ref. (34| . Near the pion threshold, the maxi- 
mum total pion angular momentum taken into account is 
up to p-wave. Therefore, all fi can be expressed in terms 
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The corresponding helicity amplitudes (Hi) are given by 

Ht = ^cos^si n ^(/ 3 +/4) , 

H 2 = -V2cos 6 -f(f x -ft- sin 2 ^(/ 3 - Si)) , 

ff 3 - -^sin^ S m^(/ 3 -/4) , 

H 4 = V2sin^(/ 1 + / 2 +cos 2 ^(/ 3 + /4)) , 

-x/Q 2 6** 
#5 = | cos^(/ 5 + / 6 ) , and 

I "-cm | ^ 

ff 6 = ^ sm ^(/ 5 -/ 6 ) . 

I ^cm | ^ 

Here, k cm is the photon momentum in the center-of-mass 
system. The structure functions that we measure can be 
expressed by these helicity amplitudes. 



1=1 



<\H 5 



\H 6 \ 2 ), 



<j tt = Re(H*H 3 - H*H^) , and 
a LT = ^ReiHKHt - H A ) + H£(H 2 



H 3 )) 



Since we are interested in the threshold region, we focus 
on the seven complex multipoles (E +, E\ + , Mi+, Mi-, 
5o+, Si+, and Si-) with I < 1, and hence 14 quantities 
must be fit to our data or determined from other data. 



In particular, the multipoles M x 



3/2 



E 3/2 



Sf^ 2 have large 



resonance contributions from the A(1232). Their Q 2 de- 
pendences at the resonance pole are well determined from 
the CLAS data and are used in this fit, whereas the 
W dependences are taken from SAID [35| . 

The other multipoles M x _, Si- and M^ 2 , e\ /2 , S\ /2 
arc fit along with the dominant Eq + , Sq+ multipoles by 
using small start values for Mi_ and Si-. Following this 
procedure, we are able to fit the measured cross section 
data and to extract the Eq + multipole. Figure |2"01 shows 
the results of the E 0+ multipole in terms of Q 2 between 
2.0 and 4.5 GeV 2 for the lowest W bin near pion thresh- 
old. 

Figure shows the Q 2 dependent Eq + multipole ex- 
traction from both our multipole fit (top) and the LCSR 
(bottom) method for three W bins with LCSR [13 and 
MAID2007 model. As soon as W is further above the 
threshold, s-wave dominance becomes weaker and reso- 
nance and higher order partial waves start to impact the 
Eq + extraction. This expected impact is clearly visible in 
our multipole fit results and is the reason why in the first 
place the measurement of the generalized form factors has 
to be carried out at the pion threshold. 
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FIG. 20: (Color online) Q 2 dependence of Eo+ normalized by 
the dipole form factor from the multipole fit. The red bullets 
are the Eo+/Gd results based on LCSR without taking the 
pion mass into account. The black bullets are the results of 
multipole. Curves as in Fig. [THJ 



VIII. SUMMARY 



We extracted the Eo+ multipole near pion threshold 
for W = 1.11 - 1.15 GeV at high Q 2 = 2.12 - 4.16 GeV 2 
with vanishing and actual pion masses and by taking into 
account different G\ form factor parametrizations. The 
results for vanishing pion mass show that Eq + /Gd is ap- 
proximately 0.3 GeV -1 and almost Q 2 -independent at 
threshold. This amplitude is larger than the MAID2007 
prediction and a little smaller than the LCSR prediction, 
which has a steeper Q 2 dependence. The Q 2 -independent 
behavior of the data may be caused by the LCSR method, 
which is based on the chiral limit — > 0. The results 
from the multipole fit method are consistent with the 
LCSR method for the lowest W bin. Independent of pion 
mass and parametrization considerations, the nn + 
channel is dominated by the transverse s-wave multipole 
Eq + . A lack of asymmetry data near the pion threshold 
does not allow us to extract the generalized form factor 
G2, but the E 0+ multipole extraction allows us to ob- 
tain Gi, and the axial form factor Ga using Eqs. (|5|) and 
(JU). Figure [22] shows the Q 2 -dependent Gi (left) and Ga 
(right) near-pion threshold. These data give strong con- 
straints on theoretical developments, especially on the 
extrapolation away from threshold and away from the 
chiral limit. 
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